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Abstract: The uncertainties involved in the calculation of bottom quark (6-quark) cross- 
sections in deep-inelastic scattering at HERA are studied in different phase space regions. 
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requiring at least one high energetic muon from the semi-leptonic 6-quark decay or a jet with 
high transverse energy are investigated. For each case the uncertainties due to the choice 
of the renormalisation and factorisation scales as well as the 6-quark mass are estimated in 
the perturbative NLO QCD calculation and furthermore uncertainties in the fragmenation 
of the 5-quark to a i?-meson and in its semi-leptonic decay are discussed. 
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1. Introduction 

Quantum Chromodynamics (QCD) has proven to give a correct description of the strong 
interaction, if a hard scale is involved in the process such that perturbation theory is ap- 
plicable. The large charm (c) and especially the bottom (6) quark mass {rrib) provide such 
a hard scale, since nif, ^ A, where A ~ 250 MeV is the typical QCD scale. The large 
mass screens the collinear singularities, so that there is no need to apply a jet algorithm 
or to subtract them into a fragmentation function in the perturbative cross section cal- 
culation. Non-perturbative contributions to heavy quark momentum distributions are of 
order 0{A/mh). The production of bottom quarks should therefore be well described by 
perturbative QCD calculations, but the inclusion of non-perturbative effects is unavoidable 
to get a physically complete result. 
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Figure 1: (a) Leading order Fcynman diagram for the reaction bbX, {b c d) some generic 

Feynman diagrams of the ©(a^) contributions. 



The 6-quark production cross-sections in strong interactions have been measured in 
proton-antiproton colhsions at SPS [|l| and at Tevatron Q-^] and, more recently, in two- 
photon interactions at LEP Q and in photon-proton colhsions (7^) at the electron-proton 
{ep) collider HERA |jlO|, |lT|. For all these measurements, except the early SPS data j^], the 
measured 6-quark production cross-sections lie above perturbative QCD expectations cal- 
culated up to next-to-leading order (NLO) in the strong coupling, a^. Since the production 
of bottom quarks is an important background process in searches for phenomena beyond 
the Standard Model at present and future high energy colliders, the understanding of the 
observed excess of the measurement over the perturbative QCD expectation, has recently 
been lively discussed. Possible proposed solutions involve higher order effects estimated 
by resummation of large logarithms [12, 13], new parton evolution schemes |14, |T6| , 
fragmentation effects [^] or even the exchange of low-mass supersymmetric particles |l8| . 
Meanwhile the last possibility seems to be excluded from LEP data [p!s|]. 

Deep-inelastic scattering (DIS) offers the unique opportunity to study the production 
mechanism of bottom quarks in a particularly clean environment where a point-like pro- 
jectile, a virtual photon (7*), collides with a proton. It will be interesting to see, if the 
excess observed in photoproduction [10, 11 [, where the photon is real and can develop a 
hadronic structure, persists for virtual photons. The large centre-of-mass energy of HERA 
colliding 27.5 GeV electrons on 920 GeV protons allows a sufficient number of b-b pairs to 
be produced and provides therefore an excellent and clean testing ground. 

The calculation of the 6-quark production cross-section requires the knowledge of the 
parton densities (mainly the one of the gluon) , of the perturbatively calculable hard parton 
parton subprocess and the correct modeling of the long range effects binding the 5-quark in 
the hadrons. The Feynman diagram of the basic production meachanism for the reaction 
7*p — > bbX is sketched in leading order (LO) of the strong coupling 0{as) in figure [l|a. As 
examples some generic Feynman diagrams of the NLO contributions are shown in figure 
—d. The dominant process is the fusion of the virtual photon emitted by the electron with 
a gluon from the proton. The production of 6-quarks is therefore directly sensitive to the 
gluon density in the proton. Furthermore the gluon evolution can be tested in the presence 
of several possible hard scales like m^, the photon virtuality or the transverse energy 
of the jet initiated by the bottom quark. 

A recent theoretical review on the NLO QCD calculation of heavy quark cross-section 
in DIS can be found in [2C[. 



In this paper we investigate the uncertainties involved in the calculation of the bot- 
tom quark production cross-section in DIS for different cross-section definitions and we 
address the question how the calculations can be best compared to experimental results: 
a short description of the perturbative NLO QCD calculation based on the HVQDIS pro- 
gram |2^, |2^, the model of the fragmentation of the 6-quark in a i3-meson and its 
semi-leptonic decay to a muon and a jet. Tools to simulate 6-quark production are de- 
scribed in section ^. After a short overview of the experimental techniques to identify 
bottom quarks, different cross-section definitions based on the measurable quantities: the 
scattered electron, the muon from the semi-leptonic fe-quark decay and the jet initiated by 
the 6-quark are presented in section |3[ Section ^ presents the necessary elements needed 
for the calculation of the cross-section for the reaction ep ^ bb ^ e njetX. The LO QCD 
calculation based on HVQDIS is compared to the LO QCD models Monte Carlo simula- 
tions in section ^. The uncertainties involved in the NLO QCD cross-section calculations 
are discussed in section |6|. This includes the uncertainties of the factorisation and the 
renormalisation scales, the 6-quark mass, the parton densities, the phenomenological de- 
scription of the fragmentation process and other non-perturbative effects. In addition, the 
uncertainties involved in extrapolating the measurable cross-section ep — > 66 — > efijetX 
to a more inclusive cross-section like e.g. ep ^ bb ^ eX are discussed in section 0. 



2. Leading order QCD Monte Carlo simulation programs 



QCD Monte Carlo models are able to simulate the hadronic final state of a DIS event 
in full detail. The four-momenta of all produced particles are made explicitly available. 
These programs are indispensable tools to correct the experimental data for detector effects. 
Although they only contain the leading order matrix elements of the hard subprocess, they 
include higher orders in a leading logarithmic approximation and they make use of models 
providing detailed treatment of the non-perturbative fragmentation phase. Here, there 
are mainly used to test the simple fragmentation model implemented in the NLO QCD 
calculation and to study the uncertainty of extrapolating measured cross-section to more 
inclusive ones. Another purpose is to study possible effects from higher orders implemented 
in the leading logarithm approximation. 

The RAPGAP Monte Carlo Q incorporates the 0{as) QCD matrix elements (ME) 
and approximates higher order parton emission using the concept of parton showers (PS) 

QCD radiation can 



1 25] based on the leading logarithm DGLAP equations |2f:, 27, |2J 
occur before and after the hard subprocess. The formation of hadrons is simulated using 
the LUND string model 1 2£ ] as implemented in JETSET pH] . As an option parton showers 



can also be simulated using ARIADNE |31]. In this model gluon emissions are treated 
by the colour dipole model |32, ^ (CDM) assuming a chain of independently radiating 



dipoles spanned by colour connected partons. 



In the CASCADE Monte Carlo simulation program |34] higher order parton emissions 

j] evolution equations are matched to a 0{as) matrix de- 



based on the CCFM H 



ment |39|-|42], where the incoming parton can be off-shell. The CCFM evolution equations 
are based on /cT^-factorisation and angular ordering which is a consequence of colour coher- 



ence, i.e. due to the interference properties of the radiated gluons. As a result in the appro- 
priate hmit they reproduce the DGLAP [|6|, 1^, |2| and the BFKL ||, |4|, |5| approxima- 
tion. At small values of the parton momentum fractions z, a random walk of the transverse 
parton momenta kx is obtained. The CCFM evolution is based on /cT^-factorisation |46|, 
where the partons entering the hard scattering matrix element are free to be off-shell, in 
contrast to the coUinear approach (DGLAP) which treats all partons entering the hard 
subprocess as massless. Off-shell matrix elements of heavy flavour lepto- and hadropro- 
duction processes have been calculated in ref. [46|. The unintegrated gluon density^ used 
in CASCADE is extracted through a fit to the proton structure function F2 measured at 
HERA [^]. The fit was performed in the range x < 10^^ and > 5GeV^. Recently, 
it has been shown that CASCADE is able to correctly reproduce the 6-quark production 
cross-section in pp-collisions at Tevatron |14]. 

In all calculations no QED corrections have been included. 



3. Measurement and definition of bottom quark cross-sections 

When a 6-quark with high transverse energy is produced in a hard strong interaction, 
a jet can be measured in the final state. The jet is composed by a hadron containing 
the 6-quark (in most cases a S-meson) and other hadrons produced in the fragmentation 
process. Due to the large 6-quark mass the jet is expected to be wider than a jet at similar 
transverse energy initiated by a light quark. However, since the rate of 6-quarks is much 
smaller than the one of light quarks, this characteristic feature is not sufficient to clearly 
identify the presence of a 6-quark . Therefore usually the semi-leptonic decay modes of 
the i?-meson to electrons and to muons B — > eX and B — > are exploited to identify 
the presence of a 6-quark. This largely suppresses the background from the light quarks 
u, d and s, which basically only can fake such a signature when a hadron produced in a 
light quark jet is misidentified as an electron or muon. The remaining background from 
charm quarks can be statistically separated from the 6-quark signal, since the lifetime of the 
meson containing the charm quark and the kinematics of the meson decay is significantly 
different. This is due to the smaller charm mass. In many cases the analysis of the muon 
channel is experimentally easier than the one of the electron channel. Therefore, here we 
only refer to the muon channel.'^ 

The experimental method often relies on the measurement of the transverse momentum 
(P^'^^) of the muon produced in the semi-leptonic decay relative to the jet axis. Due 
to the higher 6-quark mass the spectrum of this quantity is harder for jets initiated by 
6-quarks than by charm quarks. When this method is used, a jet and a muon have to be 
measured. Complementary information can be obtained from lifetime measurements of the 
-B-meson with the help of high precision silicon vertex detectors. A particularly elegant and 

Unintegrated means that the gluon density depends on the transverse parton momenta emitted along 
the cascade. This dependence is not integrated out as in the DGLAP approach where the gluon density 
only depends on the energy fraction x and on the squared transverse momentum transfer. 
^AU the arguments are, of course, also valid for semi-leptonic decays into electrons. 



simple method measures only the impact parameter of the muon track with respect to the 
main interaction vertex. Therefore, in principle in this method only a muon is needed to 
extract the 6-quark production cross-section. An alternative 6-quark tag can be obtained 
via the identification of a D*-meson. This method gives access to low 6-quark momenta, 
but is not considered here, since the efficiency of this method is rather low. 

Besides the presence of the scattered electron the experimental definition of a 6-quark 
cross-section in DIS therefore requires the presence of a muon with typically a momentum 
bigger than a few GeV and/or the presence of a jet associated to the muon at high transverse 
energy in the detector. The 6-quark cross-section measurements should always be quoted 
for the reaction ep — ^ efx{jet)X. 

To circumvent a complex cross-section definition, to facilitate the evaluation of the 
NLO QCD calculation and to make a comparison between different experimental results 
easier, in the past often extrapolations of the experimentally measured cross-section to a 
more inclusive definition have been made using leading order QCD Monte Carlo simulation 
programs. Such a procedure can only be adopted, if the models used for the extrapolation 
is reliable. Since in many cases the involved extrapolations are large, it is doubtful whether 
the accuracy of the leading order QCD Monte Carlo simulation programs is sufficient. In 
any case, unnecessary model uncertainties arc introduced in the experimental cross-section 
measurement which, when not properly stated, can make the direct comparison of data 
very difficult, if the extrapolation have been based on different model assumptions.^ 

Besides the experimentally required cross-section definition for the reaction ep ^ bb ^ 
efj,jetX (a"bjet,muon)) in tbis paper we will investigate the total inclusive 6-quark cross- 
section only defined by the scattered electron measured in the detector (cTb), the 6-quark 
cross-section where in addition a hard jet is measured ((Tbjet) and the 6-quark cross-section 
where a muon, but no hard jet is required (cJb,muon)- 

The scattered electron should be well measured in the main part of the detector. 
Therefore we require that the minimal of the exchanged photon should be > 2 GeV^ . 
Since with the present HERA data sample very high values are not accessible, is 
restricted to < 1000 GeV^. The quantity y = Q'^/{sx), where x = Q'^/{2Pq) is the 
Bjorken scaling variable, q (P) is the four- momentum of the exchanged photon (the proton 
beam) and s is the squared centrc-of-mass energy, is set to 0.05 < y < 0.7. The lower cut 
on y ensures that the squared invariant mass of the hadronic final state W'^ = (q + P)'^ = ys 
is large enough to produce a final state with high transverse energy. The higher cut on y 
is needed to ensure a good reconstruction of the event kinematics and to suppress non-DIS 
background.^ 

The jet should be well inside the acceptance of the calorimeter and the muon well 
inside the acceptance of the muon chambers and the inner tracking detectors. We require 

^The impact parameter is the distance of closest approach of the muon track to the primary event vertex. 
The average impact parameter measured in the laboratory frame is directly related to the lifetime of the 
-B-meson in its rest frame. 

*This is almost always the case, if older data are compared to recent ones. 

^High values of y correspond to low energies of the scattered electron where the unambiguous identifi- 
cation of the scattered electron becomes more difficult. 
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Table 1: Definition of the 6-quark production cross-sections used in this analysis. For each defini- 
tion the applied cut is marked with the symbol x. 

^ non-perturbative 
/ pQCD \^ fragmentation semi-leptonic 
\ function gecay 




p / \ X 
Figure 2: Sketch of the production of a J5-meson in DIS and its subsequent decay. 



here —2 < ryjet.iab < 2.5 and 30° < ^^,iab < 160° which corresponds to —1.7 < 77^,iab < 1-3. 
When the muon should be identified in the muon chambers it needs a sufficiently large 
momentum to penetrate the calorimeters. In our case we consider P^^ab > 2 GeV. 

A hard interaction is signaled by the presence of a jet with high transverse energy. The 
jet definition is based on the inclusive fcy-algorithm |4^, |4^. To cluster particles or partons 
to jets the Et scheme is used. The resulting jet four-momentum is therefore massless. This 
algorithm is well suited to make quantitative comparisons with NLO calculation, since it 
is infra-red and collinear safe |50|. To remove the purely kinematic dependence of the 
transverse jet energy on the photon virtuality Q^, the transverse jet energy is measured 
in a frame where the photon and the proton collide head on. Such a frame is e.g. the 
Breit frame defined by IxP -|- g = 0. A jet with a transverse > 6 GeV is 

required in the Breit frame. The jet should be within the detector acceptance. We require 
here: -2 < r/jet,iab < 2.5. 

The different cross-section definitions are summarised in table ||. 



4. Perturbative LO and NLO QCD calculations and fragmentation models 

The 6-quark cross-section measurements often rely on the identification of a muon from the 
semi-leptonic decay of a i5-meson. To be comparable to the measurements the calculation 
of the 6-quark cross-section proceeds in three steps as it is illustrated in figure |2|: 



1. The perturbative QCD calculation of the hard parton parton scattering cross-section 
convoluted with the parton densities and the strong coupling for the reactions ep — > 
66, ep ^ bbg and ep — > bbq 

2. The fragmentation of the 6-quark into a i?-meson described by a non-perturbative 
fragmentation function 

3. The semi-leptonic decay of the i3-meson into a muon and hadrons. 

4.1 Perturbative LO and NLO QCD cross-section calculations 

The NLO QCD heavy quark cross-sections for DIS can be calculated using the HVQDIS 
program |2^, |2^, |2^. For each event all contributions of the reactions ep bb, ep — > bbq 
and ep — > bbg are evaluated. Since the four-momenta of the outgoing particles are made 
explicitly available, the single and double differential distributions as well as correlations 
among all outgoing particles can be studied. Furthermore HVQDIS allows to easily apply 
experimental cuts. 

The 6-quark production cross-section in ep scattering can be expressed as a product of 
the strong coupling, the parton densities and the hard scattering parton cross-section which 
can be expanded in a series of perturbatively calculable coefficient functions including the 
matrix elements and the phase space factors. While the total ep 6-quark cross-section is a 
well defined physical quantity which can be compared to experimental data, the individual 
terms depend on the used calculation scheme. 

The divergences of the coefficient functions calculated to the order 0{al) due to soft 
gluon emissions are compensated by contributions from virtual gluon exchange using the 
subtraction method. The renormalisation is carried out in the modified minimal subtraction 



scheme (MS) for light quarks. The divergences from heavy quark loops are subtracted such 
that the heavy quark mass only enters in the coefficient functions. In the energy evolution 
of the strong coupling and in the evolution of the parton densities equations only light 
quarks appear. The CTEQ5F4 [^l| parameterisation of the parton densities is used as 
default. In this parameterisation the number of active quark flavours in the proton is fixed 
to four. As default, the mass of the 6-quark is set to nih = 4.75 GeV. The uncertainty on 
this value is assumed to be ±0.25 GeV. 

Uncertainties in the perturbative calculations are due to the choice of the renormal- 
isation (/x^) and the factorisation (fi'j) scales and due to the uncertainties on the 6-mass. 
In addition the uncertainty on the parton densities, in particular the gluon density, con- 
tributes. The uncertainty on the parton density function is evaluated using the ZEUS NLO 



QCD fit to recent measurements of the proton structure function |52|. 
4.2 Fragmentation of the 6-quark 

The fragmentation of the 6-quark into a S-meson is a non-perturbative process. It can be 
described by a fragmentation function obtained from comparisons with data. The frag- 
mentation of 6-quarks has recently been carefully studied in e~^e~ collisions using weakly 
decaying i3-mesons l53]-[E7i. Several fragmentation functions describing the distribution 



of the 6-quark energy fraction carried by the i?-meson have been tested and parameterisa- 
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Figure 3: Shape of the distribution of the energy fraction xb — 2EB/y/s carried by the 
_B-meson with respect to the energy of the 6-quark measured in e^e~ colhsions. Shown are data 
together with some selected fragmentation function parameterisations (see text). 

tions which describe the data best have been provided. Since in e^e~ collisions the initial 
energy of the 6-quark is known, they are usually parameterised as a function of the ratio of 
the 5-meson energy (Eb) and the beam energy (-Ebeam) which corresponds to the energy 
of the 6-quark : 

Eb '^Eb 



Xb 



(4.1) 



where ^/s is the centre-of-mass energy of the collision. The variable xb can be interpreted 
as the fraction of the 5-quark energy which is carried by the i?-meson . 

The 6-quark fragmentation functions as a function of x^ measured by the SLD, OPAL 
and ALEPH collaborations are shown in figure ^ together with some selected parame- 
39|. The data have been measured at a centre-of-mass energy which 



terisations |57, [ 

corresponds to the mass of the Z'^-boson. They are consistent with each other. The pa- 
rameterisations based on phenomenological models are roughly able to reproduce the data. 

However, the energy fraction of the 6-quark just before entering the fragmentation pro- 
cess, i.e. the 6-quark energy fraction after additional parton emissions, is not experimen- 
tally accessible (see also the sketch on top of figure ^ . The knowledge of the distribution 
of this variable, called z in the following, is, however, necessary to make predictions for 
S-meson production cross-sections in other reactions.^ 
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^In HVQDIS z is defined as pb ~ zpb, where pb (pb) is the momentum of the 6-quark (_B-meson) 
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The fragmentation function f(z) can in principle be obtained from the B-meson cross- 
section by folding them with the 6-quark production cross-section {(if,): 



as (X J abiz,npp,...) f{z,;4'p)dz, (4.2) 



where //|p is the factorisation scale which determines which contributions should be ab- 
sorbed in the perturbative part and which in the non-perturbative fragmentation process. 
The separation between the pQCD calculation and the fragmentation function is arbitrary, 
parts of the calculable contributions can be absorbed in the fragmentation function and vice 
versa. When evaluating the physically observable i?-meson cross-section the perturbative 
and the non-perturbative part have to be carefully matched to each other. In this sense 
the extracted fragmentation function always depends on the scheme how the perturbative 
cross-section ah has been calculated. As a consequence, when applying the fragmentation 
function f(z) for cross-section predictions in a different reaction, care has to be taken that 
the same approximations are applied in the perturbative calculation. In praxis, this can 
be quite tricky, since in a different reaction different parts of the perturbative calculation, 
like e.g. resummations of logarithms etc., might have to be taken into account. Since in 
e^e~ collisions at LEP or at the SLAC linear collider the centre-of-mass energy is usually 
much larger that the 6-quark mass, in the perturbative calculable part large logarithms of 
the form log (s/mf) have to resummed (NLL) |6^]. It is not clear that such contributions 
also have to be included for the cross-section calculations in pp-collisions or DIS. Here, 
rather a resummation of terms log {E?j. ^/mf^), where Et^ is the 6-quark transverse energy. 



is needed |6l| . 



As default, we follow the suggestion in refs. |jT^, ^] to use fragmentation functions 
where the fourth Mellin moment, defined as Jn = Jq z^"^ f{z)dz with = 4, folded with 
a perturbative NLL calculation is in agreement with the e"'"e~ data. Such a procedure is 
optimal when the calculations are compared to data which exhibit a cross-section behaviour 
like dab/dEx^b oc l/dE^^^, which is the case in 6-quark production at Tevatron as well 
as in DIS at HERA for large transverse momenta of the 6-quark (see line in figure ^).'' 
The parameterisation of Kartvelishvili et al. ||5^|, f{x,a) = (a + 1) (a + 2) (1 — x) is 



proposed |17, 62 1, since it is consistent with an expansion in powers of X/nib, where A ~ 0.4 



is a typical hadronic scale, as it is predicted by QCD |63, 64, 55|. The parameter^ a = 27.5 
is obtained from a fit to the ALEPH data |5^, The perturbative part has been 
calculated to NLO accuracy together with a resummation of logarithmic terms to all orders. 
It is not clear whether such a resummation is also needed for the range of Ex^b presently 
accessible in DIS at HERA. 

Therefore and also since the spectrum of the transverse 6-quark momentum deviates 
from the power-law behaviour (see above) at low values, one can use different parame- 



^The reason is that then in the moment space the 6-quark cross-section is just given as product of l/E^j, 
and the nth-moment of the fragmentation function fN{z), since 
dat/dET,bCC f dzdErfi f{z) 5{zET,b - ET,b) = fN{z). 

*Note, that this value for a is quite different then the one best describing the {l/N)dn/dxb distributions, 
namely a — 3.9 which is shown in figure ^ 



terisations with different adjustable parameters to investigate tlie uncertainty due to the 
fragmentation function. As alternative e.g. the parameterisation^ of Peterson et al. |58| 
can be used. 



4.3 The semi-leptonic decay of the S-meson into a muon and a jet 

After the i?-meson is produced it decays via a charged current interaction. In some cases 
the i?-meson decays semi-leptonically, e.g. a muon and a jet is produced. The muon can 
either be produced by a direct semi-leptonic decay {B fj,X) or by an indirect decay where 
the 6-quark decays first to a charm quark which subsequently decays semi-leptonically to 
a muon {B — > cX — > fiX). About 10% of all S-meson decays are direct decays into 
muons, in 8% (2%) of the cases the muons are indirectly produced by a charm (anti-charm) 
decay. Other decay modes of the 5-quark are much smaller, e.g. B J/tpX ^X is 
approximately 7-10~^ and B tX — > ^jlX is ~ 7-10^'^ . The branching fractions presently 
implemented in JETSET agree within 3% with the measured ones. The agreement in the 
different channels is reasonable. In the HVQDIS calculation the cross-section is determined 
using the sum of the branching ratios of direct and indirect decays of the i?-meson decays 
into a muon is = 0.22. 

In the calculation of the cross-section ep ^ bb ^ fj,X it has to be taken into account 
that the measured muon can stem from the 6-quark or the 6-quark. For this purpose, it 
is assumed that the 6-quark and the 6-quark are both independently decaying into muons 
according to the branching ratio Bfj_. We ignore the case where a 6-quark produces two 
muons via b cfiX — > fj,^X. If one of the muons is in the detector acceptance, the 
event is counted in the calculation. To improve the efficiency of the calculation, one can as 
an alternative always decay one of the 6-quarks in a muon and then decay the other one 
according to the probability that exactly two muons are in the event under the condition 
that there is at least one muon, i.e. -P2/i|>i/i = ^fi/C^^f^ ~ ^fi) — ^^J./i'^ ~ -^m)- "^^^ cross- 
section has then to be multiplied with the probability that there is at least one muon 
P>i^ = 1 — (1 — Bfj_){l — Bfj_) = 2B^ — Bf^. To ensure the cancellation of the different 
divergences appearing in the NLO QCD calculation care has to be taken that for one event 
defined, e.g. by a common x and value, only one common z-value is used for the different 
contributions appearing in a NLO QCD calculation. 

To correctly model the kinematics of the muon decay a parameterisation of the muon 
momentum spectrum in the centre-of-mass frame of the i?-meson is needed. Since the muon 
spectrum for direct (B — > fiX) and indirect {B — > cX fiX) i?-meson decays is different, 
the correct mixture has to be used. Muons are generated isotropically in the rest frame of 
the i?-meson. The absolute momentum p is obtained with the help of a probability function 
extracted from the muon momentum spectrum as obtained from JETSET. Figure |^ shows 

®The explicit form is: N^{1 — ^ — Tzr^)"^! where A'^ is a normalisation factor and e is an adjustable 
parameter, which should be in the order of e = m? /M'% ~ 0.002, where m {Mb) is a typical mass of a light 
(B) meson. Note, that the value e = 0.002 is a crude estimation and not a result of a fit to the data. 

^°The term B-meson refers here to the mixture used in JETSET, i.e. 0.4025 5° + 0.4025 B+ + 0.094 Bs + 
0.101 Ai,. This mixture is in agreement with the data (0.38 ± 0.13) B° + (0.388 ± 0.13) B+ + (0.106 ± 
0.13) + (0.118 ± 0.2) Ai,. 
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Figure 4: Momentum distribution of the decay niuons in the rest frame of tlie parent _B-meson. 
The contributions from direct (6 — > /iX) and indirect (6 — > cX — > iiX) decays are shown separately. 

the muon momentum distribution in the rest frame of the S-meson, dN/dp. The total 
distribution is shown as solid line, the distributions corresponding to the direct and indirect 
-B-decays are shown as dashed and dotted lines separately. The direct S-meson decays 
produce events with higher muon momenta. 

5. Comparison of the leading order calculations 

To check that the simple form of 6-quark fragmentation and of the subsequent decay of the 
i3-meson implemented in HVQDIS gives reasonable results, a comparison to the more com- 
plex fragmentation models used in the Monte Carlo simulation programs can be made. For 
this purpose, the same free parameters in the perturbative cross-section calculations have 
to be chosen. A comparison of the LO QCD prediction of the inclusive DIS 6-quark cross- 
section calculated by HVQDIS and by RAPGAP^^ is shown in figure ^. The renormalisation 
and factorisation scale = /i^ = //^ = has been used and the 6-quark mass was set 
to rrih = 5 GeV. Moreover the 1-loop formula of the strong coupling Og and the CTEQ5L 
parton density functions have been used. 

^^To obtain a pure LO calculation in RAPGAP, the parton shower option has been switched off. Moreover, 
for the calculation of the differential 6-quark and jet cross-section the fragmentation has been switched off. 
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Figure 5: Inclusive DIS &-quark cross-sections as a function of (a) the transverse energy of the 
6-quark -ET,b,Brcit, (b) the transverse energy of the jet -ETjot.Brcit, (c) the transverse energy of the 
muon -ET,/i,iab, (e) the pseudo-rapidity of the jet Tyjot.iab, (d) the pseudo-rapidity of the 6-quark 77b, lab, 
(/) the pseudo-rapidity of the muon ry^^jab- Shown are LO QCD cross-section calculated by HVQDIS 
(LO) and by RAPGAP without parton showers. More details are given in the text. 

The inclusive DIS 6-quark cross-section as calculated by LO-HVQDIS is Ub = 530 pb 
and agrees within 2% with the prediction from RAPGAP. For the 6-quark cross-section 
requiring a jet dbjet = 148 pb is found. The agreement with RAPGAP is within 4%. 
This is a strong indication that the LO matrix elements used in the two programs are 
the same. However, difference up to 5% are found for large transverse energies of the 
6-quark i?T,b,Breit and of the jet -Er,jet,Brcit- This can be seen in figure |5| where the inclusive 
6-quark cross-section as a function of the transverse energy of the 6-quark ^'r.b.Breit (a) 
and as a function of the transverse jet energy -E^T',jet,Breit (t*) the Breit frame is shown. 
No strong dependence on the pseudo-rapidities of the 6-quark ?7b,iab or of the jet pseudo- 
rapidities in the laboratory frame ?7jet,iab is found (see figures |d and |e). 

For the 6-quark cross-section requiring a jet and a muon a"bjet,muon = 19 pb is found 
with LO-HVQDIS, while the one from RAPGAP is 10% larger. The observed difference 
in (Tbjet.muon is probably due to an incomplete modeling of the muon fragmentation in the 



HVQDIS program. This can be seen in figures and |5[/. The transverse energy spectrum 
of the muon is harder for the HVQDIS calculation. At large transverse energies it is 10% 
higher. The biggest difference is found in the pseudo-rapidity spectrum. In the RAPGAP 
calculation the muon is produced more forward. Difference up to 30% are found. A possible 
explanation of this effect is that the fragmentation model in HVQDIS is too simple to cope 
with the multi-parton environment where it is important that all colour connections are 
correctly defined. 

6. Perturbative NLO QCD calculations and their uncertainties 
6.1 Calculated NLO QCD cross-sections 

The NLO QCD 6-quark cross-sections for the four cases a^, o"b,jeti fb.muom crb,jet,muon are 
given m table |. The cross-sections have been calculated with the renormalisation and 
factorisations scale set to //^ = ^ -|- 4m^, where pT,b is the transverse momentum of the 
6-quark in the Breit or hadronic centre of mass frame. The 6-quark mass has been set to 
nib = 4.75 GeV. The CTEQ5F4 parameterisation has been used for the parton densities. 
Also shown are the charm quark cross-sections.^'^ 

The inclusive charm quark cross-section is about 370 times larger than the inclusive 
6-quark cross-section. This is also expected, since the electrical charge of the charm quark 
is higher and its mass is smaller. The heavy quark mass enters directly in the matrix 
element and in the phase space factor, in addition it changes the parton kinematics such 
that different parton densities are probed. It is, however, interesting that if a hard jet or a 
muon is required in the detector acceptance the charm quark cross-section is only an order 
of magnitude bigger than the 5-quark cross-section. The reason for this behaviour is that 
the large 6-quark mass naturally produces particles at higher transverse energy. 

The measurable 6-quark cross-section (crb,jet,muon) is about 20 times smaller than the 
inclusive one (fib). Only part of this difference is due to the branching fraction. The 
large fraction of the 6-quark cross-section can not be measured in the detector, since the 
6-quark has either a low transverse energy or is produced in the forward direction^'* outside 
the detector acceptance. 



NLO 


cJb (pb) 


c^bjet (pb) 


<7b,muon (pb) 


'^bjet.muon (pb) 


6-quark 


598 


230 


53 


26 


c-quark 


24170 


2320 


515 


188 



Table 2: NLO QCD cross-sections calculated for the four studied cross-section definitions for 
6-quarks and c-quarks. 

^^The hadronic centre of mass frame is defined as g* -I- P = 0, where q (P) is the photon (proton) 
momentum. The centre of mass frame can be transformed by a longitudinal boost to the Breit frame. The 
transverse energy is therefore the same in both frames. 

^^The charm quark cross-section have been calculated with the factorisation and renormalisation scales 
set to = (. -I- 4m^ and the charm quark mass is set to rric — 1.4 GeV. The CTEQ5F3 parton density 
functions and the Peterson fragmentation function with e — 0.033 have been used. The semi-leptonic muon 
momentum distribution and the branching fraction Bc^^x have been modified for charm quarks. 

^*At HERA, the proton moves into the -I- z- direct ion. The forward direction is the therefore the region 
located toward the proton remnant. 




Figure 6: Cross-sections (ub, crbjot, fb.muon, o'bjct,muon) as a function of the transverse energy (a) 
and the pseudo-rapidity of the 6-quark for the four studied definition (see text). The inner band 
represents the renormalisation and factorisation scale uncertainties and the outer band in addition 
the uncertainty due to the 6-quark mass. The renormahsation and factorisation scale has been set 
to /^^ = J, -|- Ami. ^ comparison, the solid line indicates the power law ab/dEx.b oc E^^. 

The 5-quark cross-section as a function of the 6-quark transverse energy in the Breit 
frame and its pseudo-rapidity in the laboratory frame is shown in figure ^. The central 
NLO QCD predictions are shown as lines. The inner band indicates the renormalisation 
and factorisation scale uncertainties, the outer bands indicates in addition the uncertainty 
due to the 6-quark mass. For the inclusive 6-quark cross-section (cTb) the 6-quark transverse 
energy distribution exhibits a broad peak at about Ex^h^Breit ~ 

5 GeV. Above this value 

the distribution falls like da\^/dExfi^BreU -S'Tfe Breit toward larger values. Below this value 
the 6-quark cross-section is approximately constant. If a hard jet is required (dashed or 
dashed-dotted line in figure the behaviour at large ET,b,Breit is not modified, but the 
region of ET,b,Breit is suppressed. When measuring the 6-quark cross-section via the semi- 
leptonic muon decay, the NLO QCD predictions for relatively large Ex^b,Breit are tested. 
For all shown cross-sections the uncertainties at low transverse 6-quark energies are larger. 
Figure p> shows the 6-quark pseudo-rapidity in the laboratory frame. For all cross-section 
definitions most of the 6-quarks are produced centrally, i.e. —2 < %^;ab < 2. However, in 
the inclusive case (ub) about 10% of the 6-quarks are produced at pseudo-rapdities beyond 
flh,iab > 2.5, where there is no detector coverage. If the muon is required to be in the 
detector acceptance, only 6-quarks within — 2 < "qb^iab < 2 contribute to the cross-section. 

6.2 Factorisation and renormalisation scale and &-quark mass uncertainties 

A fixed order QCD cross-section calculation usually depends on the renormalisation and 
the factorisation scales. The renormalisation procedure removes divergences due to loop 
contributions in which virtual exchanged partons can have very large momenta (ultra-violet 
divergences) . The factorisation scale separates the perturbatively calculable short distance 
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Figure 7: (a) Ratio of the cross-sections (Tb, Ub.jot, crb.muon, Cbjct.muon calculated with the renor- 
malisation and factorisation scales set to /i^ = ^{Q^ + Ami) ^^'^ A*^ ^ Co(Q^ + 4m^) where = 1- 
(&) Ratio of the cross-sections calculated at NLO and at LO (K-factor) as a function of the scale 
factor ^. 



contributions from the non-perturbative long-distance contributions. The non-perturbative 
part is absorbed in the parton density functions. In this way divergences due to collinear 
parton radiation are removed. The choice of the renormalisation and factorisation scales 
in the calculation is arbitrary, but should correspond to the hard scale (Q^) involved in the 
hard scattering process. In a complete calculation a physical observable does not depend 
on the choice of the renormalisation and factorisation scale. The residual scale dependence 
in a fixed order calculation can be used to estimate the uncertainties due to neglected 
contributions. 

In a LO calculation the calculated cross-section depends on the renormalisation scale 
via the energy behaviour of the strong coupling as{^r)- When the renormalisation scale 
is increased, it is therefore expected that the cross-section decreases in the same way. In 
a NLO calculation the dependence of can be compensated by the perturbatively 
calculable coefficient functions, since they also depend on /i,.. The scale dependence can 
therefore be significantly reduced when the NLO corrections are included in the calcula- 
tion. In some cases the residual scale dependence is closely related to the higher order 
contributions which have not been calculated. 

The dependence on the scale factor ^, defined as = ■CQ^, where Q is the hard scale, is 
shown in figure for the four different cross-section definitions as calculated in NLO QCD. 
The hard scale has been chosen to be ^ = Q^-|-4m^. The cross-sections are normalised to 
the cross-sections obtained for C = Co = !> i-e. the ratio {da / d£^) / {da / d^Q) is shown. When 
^ is varied by about a factor of 10 in both directions, the inclusive cross-section Cb varies 
by about ±5%. Once a hard jet is required, the cross-section variation increases to ±20%. 



Figure |^ shows the NLO correction as a function of the scale factor ^. The NLO 
correction is defined as the ratio of the cross-section calculated at NLO to the one calculated 
at LO. The dependence of the NLO correction on the scale factor is rather small for all 
studied cross-section definitions. While the NLO correction is only about 10% for the 
inclusive cross-section cjb and for the one requiring a muon (Tb,muon) it increases to almost 
a factor of 2 when a jet is required. 

This together with the stronger scale dependence is an indication that for the cross- 
section measured via the P^"^^ -method where a jet has to be required higher order contri- 
butions are more important. However, whenever a hard jet is involved in the cross-section 
measurement, the increased uncertainty is unavoidable. It is interesting that the require- 
ment of a muon in the detector acceptance, does not lead to an increased uncertainty. This 
is clearly an advantage for cross-section measurements based on the i?-meson lifetime using 
displaced vertices. However, also in this case the requirement of a hard jet in the Breit 
frame might be necessary in order to efficiently reject the large background from inclusive 
DIS events. 

Here and in the following the renormalisation and factorisation scales have been set 
equal. Since this is an assumption with no deeper justification we have, for = ^+4m^, 
also varied /i^ and //j independently. For all factorisation (renormalisation) scales the 
cross-section rises for decreasing (increasing) //^ (M/)- '^^^ 5-quark cross-section changes 
from Ub = 597.6 pb for n'j = ^ir = Pt b ~^ = 671.3 pb for fi'j = ^{p\^ + Ami) 

and /u^ = 0.25 (p^,^ + 4m^) and to a\, = 511 pb for = 0.25 (p|,^ + 4m^) and = 
4 (p^ ^ + 4m^). Similar conclusion hold for the other cross-section definitions. We find that 
the variation with the renormalisation scale is about equal for all choices of the factorisation 
scale, i.e. about ±5% for Cb and about ±20% for (Tb,jet,muon- Also the dependence of the 
factorisation scale is about equal for all choices of the renormalisation scale, i.e. about 
±5% for (Tb and also for cibjet.muon- The fact that the renormalisation scale dependence 
does not depend on the factorisation scale choice and vice versa, remains valid even for 
extreme scale factors like ^ = 40. In the //j — ^1 plane we therefore have not found a 
plateau region where the scale dependence is significantly reduced with respect to other 
regions. 

Figure |8| shows the cross-sections fib, o"b,jet5 crb,muon) Cbjet.muon as a function of 
and of ET,jet,Breit- From 2 < < 1000 GeV^ the cross-sections fall by about 4 or- 
ders of magnitudes toward larger (figure ^). The dependence is similar for all 
studied cross-section definitions. At low Q^, ^rbjet (c'bjet.muon) is about a factor 3 (2) 
smaller than the inclusive cross-section (o"b,muon)i at ^ 200 GeV^ the ratio o"b/o"bjet 
(cb,muon/<7bjet,muon) is Only 2 (1). The dependence on the transverse jet energy is simi- 
lar for all cross-section definitions for ET,jet,Breit > 6GeV (see figure ^). The inclusive 
cross-section as well as the one requiring a muon (Xb^muon result in a mean transverse 
jet energy of about 10 GeV. When a hard jet is required the mean transverse jet energy in- 
creases to about 15 GeV. As a consequence possible terms of the form log {E"^ -^^ Breitl'^'b) 
increase from about 0.5 to 1. About 35% of the inclusive 5-quark events have a hard jet 
with Et Jet, Breit > 6 GeV, when a muon in the detector acceptance is required 45% of the 
events have a hard jet. 
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Figure 8: Cross-sections (ffb, fbjet) Cb,muon> Cb,jet,muon) as a function of (a) and the trans- 
verse jet energy -ETjot.Broitl^)- The inner band represents the renornialisation and factorisation 
scale uncertainties and the outer band in addition the uncertainty due to the 6-quark mass. The 
renormaUsation scale has been set to fi^ — 



Ami. 



The uncertainty on the renornialisation and factorisation scales can be estimated by 
varying the scale factor ^. It is common practice to vary /x^ by approximately one order 
of magnitude, i.e. the scale factor is varied by 0.25 < ^ < 4. In addition, the uncertainty 
introduced by the limited knowledge on the 6-quark mass can be estimated by varying it 
between 4.5 < nii, < 5.0 GeV. The default value is set to mi, = 4.75 GeV. 

Figure ^ shows the measurable 6-quark cross-section crbjet,muon as a function of Q^. 
The hard scale is set to = ^ -|- 4m^ . The inner band illustrates the scale uncertainty, 
the outer band shows in addition the uncertainty due to the 6-quark mass. The size of the 
uncertainty is better visible in figure ^6 where the relative difference of the modified cross- 
section calculation to the default cross-section calculation, i.e. (crbjet,muon(Co i ^,mfe ± 
Amfe) - cTb,jet,muon(Co,"2-6))/o-bjet,muon(?o,"^6), IS shown. The total Uncertainty is about 
±15% at low and slightly decreases towards larger Q^. 

Also shown in figure |^ is the pure LO calculation (dashed line) and the calculation 
complemented with parton showers (dotted line). When parton showers are included in 
the calculation, the LO cross-section increases by about 15% at low and by about 10% 
at large Q^. At low Q^, the NLO correction increases o"b,jet,muon by almost a factor of 2, 
toward larger the NLO correction increases CTbjet.muon by 1.6. When parton showers are 
included the NLO correction is lower, <Tbjet,muon increases by 1.75 at low and by 1.5 at 
large Q^. 

For 6-quark production in DIS there are three possible choices for the hard scale Q^: 
m^, and the transverse momentum of the 6-quark or of the associated jet in the Breit 
frame. The transverse momentum can either be defined to be the transverse momentum 
of the 6-quark px^b or the maximal transverse momentum of the 6-quark or the 6-quark 
PT,b,max- In principle, any function of these scales is a possible choice. If one of the possible 
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Figure 9: (a) Bottom quark cross-section fTbjet.muon as a function of for different choices of the 
renormalisation and factorisation scales /i^. (6) Scale uncertainty for the choice = Pt6 + 4to^ 
(band) and the ratio of the cross-sections when calculated with different scale choice normalised to 
the one calculated with ^ Pt b ^ Anr^ (lines) . The inner band represents the renormalisation 
and factorisation scales uncertainties and the outer band in addition the uncertainty due to the 
6-quark mass. 



scales is much larger or smaller than the others, but still sizable, it can be expected that 
large logarithms of the form log (Qf /Q^) appear in the calculation. These logarithms have 
to be resummed to make the calculation reliable. 

The relative difference of the cross-section calculated with different scale choices, i.e. 
(crb,jet,muon(Q?) " o-b,jet,muon( Qo))/^b,jet,muon (Qq) ' is shown in figure |6 as a function of 
The default scale choice is Qg = b + ^m^. Shown as examples are Qf = + 4:ml (solid 
line), Ql = Ami (dashed line), Qf = p^^^ + + Ami (dotted line), Qf = PT,b,ms.x + ^"^6 
(dashed-dotted line). All studied possible scale choices lead to cross-section predictions 
within the uncertainties of the default choice. The choices Qf = f, + + 4m^ and 
Ql = Pxb max+^"^b §1^6 exactly the same result. The cross-section calculated for = Ami 
is 10% higher over the full Q"^ region, the one calculated for Qf = + Ami higher 
at lower Q^, but only 2% higher at large Q^. 

The NLO QCD cross-section prediction is therefore more or less independent on the 
choice of the hard scale. One way to decide which scale choice should be used as the default, 
is to study the uncertainties introduced by the various possible scale choice. Figure |lO| shows 
the cross-section uncertainty for the six studied scale choices as a function of Q^. For all 
choices the total uncertainties stay within ±20%. The uncertainty is larger at low for 
all choices and decreases toward larger Q^. When the transverse 6-quark momentum is 
included in the definition of Q^, the uncertainty is in general smaller. For instance, for 
Qf = + Ami uncertainty is around ±(10 - 20)%, while for Qf = p'^^ + Q'^ -\- Ami 
uncertainty is around ±(10 — 15)%. For Ql = Ami uncertainty is ±(15 — 20)%, while 
for Ql = p^^ + Aml the uncertainty is only ±(10 — 15)%. No big difference is seen between 
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Figure 10: Uncertainty of the 6-quark cross-section o'b.jct,inuon as a function of for different 
choices of the renormalisation and factorisation scales. The inner band represents the renormalisa- 
tion and factorisation scale uncertainties and the outer band in addition the uncertainty due to the 
6-quark mass. The scales are explained in the text. 



choices involving pT^b or PT,b,i 



The smallest scale uncertainty is obtained for 

^ PT,b 



PT,b + 



Ami Qi ~ PTb max ^iT^l' "^^^ choice Qi = Pxb^ ^"^6 therefore proposed as default. 



Figure shows the inclusive 6-quark cross-section (left) as a function of the trans- 
verse energy of the jet with the highest transverse energy {ETjet,Breit) in the event and 
as a function of the pseudo-rapidity of the jet which contains the muon from the semi- 



leptonic &-quark decay. The scale is set to fj? = p^b^ 4m^. In the right part of figure 11 



the uncertainty of the NLO calculation is explicitly shown. It is interesting that the scale 
dependence of the cross-section varies as a function of the shown observables. It is most 
pronounced toward larger transverse jet energies. This is the case where the transverse 
jet energy and consequently the transverse 6-quark momentum are very different from the 
6-quark mass. A large scale dependence is also seen in the forward direction. For low 
where ET,jet,Breit ~ iT^b, the scale dependence is rather small (down to ±5%). 
This is the region where most of the inclusive 6-quark events are produced. Therefore the 
small scale dependence of the inclusive cross-section cTb which was discussed in the context 
of figure 1^ seems to be purely accidental. 
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Figure 11: The cross-section a^, as a function of tlie jet witli tlie liigliest transverse energy 
£'T,jct,Bicit(a) and the pseudo-rapidity ?7muonjct,iab of the jet associated with the muon (c). The 
inner band represents the renormalisation and factorisation scale uncertainties and the outer band 
in addition the uncertainty due to the 6-quark mass. The uncertainty, i.e. the ratio of the cross- 
section calculated with the default settings over the one with varied parameters is shown in (6) and 
(d). The renormalisation and factorisation scale have been set to /i^ = + 47ti^. 
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The total uncertainty including the 6-quark mass is about constant namely around 
±20%. Only in the large ET,jet,Breit region the uncertainty is increased up to about 
30%. In the region where the scale total uncertainty is small, the uncertainty due to 
the 6-quark mass is large. 

The increase of the uncertainties toward low ETjet,Breit is most probably due to an 
infra-red sensitive phase space region In fact, the 6-quark cross-section steeply falls 

toward ETjet,Breit ^ and even gets slightly negative in the extreme limit (not shown). 
This reflects an incomplete cancellation of the (positive) divergences of the real infra-red 
parton emissions and the (negative) divergences of the virtual contributions. The only pos- 
sibility that no hard jet is produced in the event, is the configuration where the 5-quark and 
the 6-quark are scattered exactly in (or anti-parallel to) the direction of the proton. 

Note, that the 20% total uncertainty on the inclusive cross-section fib is larger than 
the one on the cr^ jet, muon 

cross-section which is overall only about 15% for /i^ = Ptb~^ ^''^1 
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Figure 12: Same as figure |ll|, but for the measurable cross-section ah. 



jct.muon ■ 



(see figure ■ The differential cross-section for the case crb,jet,muon is shown in figure 12 
The uncertainties exhibit a similar behaviour than in the inclusive case, i.e. the scale 
dependence gets larger toward larger transverse jet energies and toward the forward region, 
while the 6-quark mass uncertainty is largest at low transverse jet energies and in the 
backward region. However, the overall uncertainty is lower than in the inclusive case. At 
low Exjet.Breit, it is Only about ±(10 — 15)% and increases to about ±20% toward larger 
ET,jet,Breit- The Uncertainties in the backward and forward region toward the end of the 
detector acceptance are about the same in the case of db and o"bjet,muon- 

6.3 Parton density uncertainties 

The calculation of cross-sections for observable processes generally involve the convolution 
of the perturbatively calculable coefficient functions corresponding to the hard subprocess 
with the non-perturbative parton densities which describe the structure of the proton. 
The parton densities absorb the collinear divergences of initial state parton radiation ap- 
pearing in the perturbative calculation. The boundary defining which contributions are 
perturbatively calculated and which are absorbed in the parton densities is given by the 
factorisation scale. The parton densities can be extracted in a fit procedure by compar- 
ing cross-sections calculated by convolution them with the coefficient functions to data. 
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Figure 13: Uncertainty on the CTbjot.muon cross-section introduced by the uncertainty on the gluon 
density as a function of (a) and -ETjot.Brcit- The inner band represents the renormahsation 
and factorisation scale uncertainties and the outer band in addition the uncertainty due to the 
&-quark mass. Shown as line is the ratio of the central ZEUS NLO QCD fit and the CTEQ5F4 
parameterisation. 

Usually, a variety of data such as inclusive DIS, Drell-Yan, prompt photon, VF-boson and 
inclusive jet production in ep or pp collisions etc. are used to optimally constrain the par- 
ton densities. Recently, significant progress has been made in determining in addition also 
the uncertainties on the parton densities. This is achieved by a careful error propagation 
in the fit procedure. 

To estimate the impact of the limited knowledge of the proton parton densities on 
the calculation of the 6-quark cross-section, we use the recent NLO QCD fit performed by 
the ZEUS collaboration Since the large part of 6-quark event are initiated by gluons, 
the cross-section uncertainty is mainly due to the gluon density for gluon energy fractions 
Xg in the range^^ 10~^ — 10~^. The hard scale where the parton densities are probed lies 
typically within 10 - 100 GeV^. 

The uncertainty of the measurable cross-section (Tbjet.n 



is shown function of 



Q in figure 13a and as a function of ETjet,Breit in figure 13b. The band indicates the 



cross-section uncertainty due to the parton densities. The solid line gives the ratio of the 
cross-section calculated with the central ZEUS parton densities and the CTEQF4 parton 
densities. The cross-section calculated with CTEQ parameterisation is overall about 10% 
lower. The smaller cross-section is only observed at low ETjet,Breit, where the bulk of 
the events are, for high ETjet,Breit the two different parton density parameterisations give 
about the same results. The cross-section difference is not bigger than the uncertainty 
determined for the ZEUS parameterisation which is not bigger than ±(5 — 10)%. 



About 90% of all 6-quark events have energy fractions Xg in this range. Due to kinematic reasons, 
values Xg < 10~^ can not be reached. About 10% of the events have Xg > 10^^. In some case also very 
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6.4 Fragmentation function uncertainties 



The uncertainty due to the fragmentation function is estimated by changing possible choices 
of the fragmentation functions and their free parameters. Since care has to be taken that 
the fragmentation function is well matched to the perturbative calculation (see section ^) , 
the estimated uncertainty is probably bigger than the real one. 



Figure 14 shows the transverse energy and pseudo-rapidity distribution of the muon 



for the inclusive case. Shown as band is the NLO QCD cross-section crt calculated with the 



Peterson fragmentation function |55] with different parameters, i.e. 0.001 < e < 0.005. The 



dashed line indicates the cross-section obtained with the parameterisation of Kartvelishvili 



et al. [£9|| using the SLD fit parameter and the dotted line is obtained by using the fit 



parameter parameter proposed in refs. |17, 62]. While for low transverse muon energies all 



cross-section agree, noticeable difference are found for ET,muon ^ 2 GeV. No big differences 
are found in the pseudo-rapidity distribution of the muon. 

Since most of the 6-quark events produce muons with low transverse energies, the total 
measurable cross-section cr^ jet, muon varies only by about 10% when different 6 parameters 
are used for the Peterson fragmentation functions and by 5% when different fragmentation 
functions are used. 

When the 6-quark cross-section db is calculated in leading order, the ET,muon spectrum 
agrees with the one calculated with RAPGAP when e = 0.009 is used. The pseudo-rapidity 
distribution, however, always peaks more forward. 




Figure 14: The 6-quark cross-section Ub as a function of the muon transverse energy (a) and 
pseudo-rapidity (b) calculated with different fragmentation functions. Shown as band is ah calcu- 
lated with the Peterson fragmentation function for e — 0.001 — 0.005 and the parameterisations 
of Kartvelishvili et al. |Q using the SLD fit parameter and the one obtained by a fit to the fourth 
moment of the 5-quark fragmentation function measured in e~^e~ collisions |l^, p2]. 



7. QCD model uncertainties 



7.1 Extrapolations uncertainties due to QCD models 

QCD models calculating the leading order 6-quark cross-section using the leading order 
matrix elements and models to implement higher order radiation are the only calculations 
presently available that provide the full hadronic final state of an event. Therefore they 
are the only available tools to study the detector response and to determine the detector 
acceptance. Often they are also used to extrapolate from the measured cross-section to a 
more inclusive cross-section. 

In table ^ the cross-section obtained using the QCD models is compared to the LO and 
NLO calculation. ME+PS denotes the parton shower model, CDM the colour dipole model 
and ME*+CCFM combines the off-shell matrix elements with initial parton radiation based 
on the CCFM evolution equations. The physics content of these QCD models is explained 
in more detail in section 

For the inclusive 6-quark cross-section cjb the same result is obtained for the ME-I-PS 
and for the CDM model. Since the inclusive cross-section does not depend on the treat- 
ment of the hadronic final state this is also expected. It it, however, interesting that the 
6-quark cross-section requiring a jet or a muon in the final state, are also very similar for 
both models. The inclusive cross-section Ub is 25% lower than the NLO QCD result. This 
difference between the QCD model and the NLO QCD calculation increases to about 40% 
when a jet or a muon is required in the cross-section definition. While for Cb the QCD 
models agree with the LO QCD calculation, they lead to a 20% — 40% bigger cross-section, 
when a jet is required. The increase is due to the inclusion of the parton showers. 

Since in the ME*+CCFM model a different gluon density in addition to the dif- 
ferent parton radiation pattern and the off-shell matrix elements is used, the inclusive 
6-quark calculated cross-section can be different. The ME*+CCFM model gives a cross- 
section which is 2 — 3 times higher than the one calculated with the ME+PS or CDM 
models and about 20% - 40% higher than the NLO QCD result. 

When the QCD models are used to remove the jet in the cross-section definition, i.e 



Model 


(Tb (pb) 


cjbjet (pb) 


0"b,muon (pb) 


fbjet.muon (pb) 


NLO 


597.6 


230.4 


52.7 


26.1 


LO 


449.4 


113.6 


40.8 


14.4 


ME+PS 


475.2 


164.7 


38.7 


17.1 


ME only 


475.2 


123.3 


38.5 


13.8 


CDM 


475.2 


163.5 


37.0 


16.3 


ME*+CCFM 


822.0 


364.4 


65.8 


35.5 



Table 3: Calculated fe-quark cross-sections for the four studied definitions. The calculations are 
made with /j? = p|, ^ + 4to^. See the text for more explanations. 



quoting o"b,muon instead of o"bjet,muon) the following extrapolation factor has to be calculated: 



MC 

data data _ b 

"-^b b,jet,muon ' ]VIC 

'^b,jet,muon 

For MEPS and CDM this extrapolation factor is about 2.3. However, for NLO it is only 
2.0 and for ME*+CCFM it is 1.8. Since it is not known, which of the calculation gives the 
correct result, the jet extrapolation introduces a model uncertainty of about 20%. If one, 
for instance, assumed that the experimental measurement would give o'bjet,muon= 35.5 pb 
as predicted by CASCADE and one used the MEPS Monte Carlo simulation to quote 
a cross-section (Tb,muon) the result would be o"b,muon = 78.1 pb. This is, however, about 
20% higher than the correct result. Moreover, the extrapolation factor is dependent on 
the choice of the renormalisation and factorisation scales. This is illustrated for the NLO 
QCD calculation in figure ^. The ratio of the dashed-dotted and the dotted line is the 
extrapolation factor needed to correct for the jet acceptance. The extrapolation factor 
varies from 1.8 obtained for scale choice of ^ = 0.1 to 2.25 for ^ = 10, which corresponds 
to a change by 20%. 

Sometimes the jet definition is kept in the cross-section definition, but the requirement 
of a muon is removed in order to get a cross-section definition where no fragmentation model 
is involved in the NLO QCD cross-section calculation. The idea of such a procedure is that 
the complex fragmentation model as implemented in the QCD models are more reliable 
than the very simple one used in the NLO QCD calculation.^^ If the muon requirement 
is removed from the cross-section definition, the following extrapolation factor has to be 
calculated: 

MC 

data _ data _ b.jet 

h,}et "bjet.muon ' MC 

b,jet,muon 

The correction factor calculated with MEPS, CDM or ME*+CCFM is about 10. The one 
calculated with NLO QCD is 8.8. The introduced model uncertainty is about 10%. The 
correction factor is rather large, only 10% of the quoted cross-section result is actually 
measured in the detector. However, part of the correction is due to the muon branching 
fraction which is well known. About 75% of the muon in the rapidity range measured by 
the detector —1.7 < r?muon < 1-3 and about 30% of the muons satisfy Pmuon > 2GeV. 

The model uncertainties of the extrapolation factors are related to the different distri- 
bution of the variables defining the cross-section, i.e. the muon and jet transverse energies 
and pseudo-rapidities. The shapes of these distributions for inclusive b-quaik events are 
shown in figure |l^. Figure [T5| a shows the transverse energy of the jet with the highest Et 
in the Breit frame and figure ^5^ shows the pseudo-rapidity of the jet which is associated 
with the muon.^^ Figure 15c and figure show the transverse energy and the pseudo- 



rapidity of the muon. Given the very different treatments of the underlying partonic final 



^^Even if this idea sounds reasonable one has to be careful, since the different partonic final state can 
lead to different distributions of the variables defining the muon kinematics. 
^'^The association of the muon and the jet is performed by the jet algorithm. 
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Figure 15: Comparison of the shape of the transverse energy and the pseudo-rapidity distribution 
of the jet and the muon for the inclusive 6-quark cross-section as obtained by the NLO QCD 
calculation and the QCD Monte Carlo models. 

state it is remarkable that the jet and muon distribution of the tested calculations are not 
too different. 

The CDM and NLO QCD have a similar transverse jet energy spectrum. The MEPS 
model produces less events at high Et- The spectrum calculated by ME*-|-CCFM is sig- 
nificantly harder. Less events at low Et and more events at large Et are expected. The 
pseudo-rapidity spectrum of the jet associated to the muon is similar for the MEPS and 
the CDM model. NLO QCD has a spectrum which is clearly shifted toward the forward 
direction. In between these two models lies the ME*-|-CCFM prediction. In contrast to the 
transverse jet energy NLO QCD exhibits the hardest transverse muon energy spectrum. 
MEPS and CDM have similar spectra. ME*-|-CCFM lies in between. At low transverse 
muon energies all calculations more or less agree. The muon pseudo-rapidity spectrum 
is most shifted toward the forward region for the ME*-|-CCFM prediction. NLO QCD is 
most backward and CDM agrees more or less with MEPS. 

The detailed understanding of these distributions is rather difficult. While the jet dis- 
tributions arc a reflection of the different partonic flnal state, the muon distributions also 
include hadronisation effects. It is e.g. remarkable that the muon pseudo-rapidity distribu- 
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tion is most shifted toward the backward region for the NLO QCD calculation, although 
the jet pseudo-rapidity distribution is shifted most forward. This might be explained by 
the observation made in the context of figure ^ that the hadronisation pulls the muon 
forward. 



7.2 Hadronisation corrections uncertainties 

If a jet is required in the cross-section calculation, the NLO QCD prediction should be 
corrected for hadronisation effects, since the jets can only be formed by the three partons 
available in the final state. Such a correction is, however, only meaningful, if the shape of 
the transverse energy and pseudo-rapidity distributions of the parton jets calculated with 
the QCD models is not too different from the one of the NLO QCD calculation. 

The cross-section fTb,jet,muon decreases by about 10% at low and by about 5% 
at high Q^, if hadronisation effects are included. The dependence of the hadronisation 
corrections on the transverse jet energy and the jet pseudo-rapidity for the cross-section 



Ubjet.muon IS shown in figure |16| a and figure |16| c. The solid line is the MEPS expectation 



based on the hadronic final state, the dashed line is the one based on the partonic final 
state. Hadronisation effects lower the cross-section at low Et, but increase the cross-section 




Figure 16: Transverse jet energy and pseudo-rapidity for the measurable cross-section 
fb,jct,muon based on the hadronic (cThad) and the partonic (upar) final state. 



at high Et- The hadronisation correction factor, i.e. cihad/'^par is shown in figure |l6p and 
figure 16 i for the three discussed QCD models. The symbol cJhad denotes the 6-quark cross- 
section calculated for hadrons, fipar denotes the 6-quark cross-section calculated for partons. 
The biggest model uncertainties on the hadronisation correction is seen at large transverse 
jet energies and for low jet pseudo-rapidities. This are, however, the regions where the 
cross-section is rather small. 



8. Conclusions 



The uncertainties involved in the NLO QCD calculation of 6-quark cross-sections have been 
estimated. Besides the inclusive 6-quark cross-section, cross-section definition requiring the 
muon from the semi-leptonic 6-quark decay or the jet induced by the 6-quark have been 
studied. The uncertainties due to the renormalisation and factorisation scale are about 
10% — 20%. The uncertainties introduced by possible scale choices are within this mar- 
gin. In general, the total uncertainties including the uncertainties on the 6-quark mass are 
about constant over the full jet transverse energy and pseudo-rapidity range. The total 
uncertainties are smallest for the cross-section where the muon and the jet are within the 
detector acceptance. When extrapolating the measurable cross-section to more inclusive 
cross-section definitions large extrapolation factors can be involved. The model uncertain- 
ties on these factors are estimated to be 20%. Hadronisation corrections lead to model 
uncertainties of about 10%. 
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